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I. INTRODUCTION 



Quantum effects in accelerators have been of interest for many years, [1] , [2] , 
and [3] ; however recently with the development of multi TeV colliders and interest 
in large linear colliders, they have become the subject of wide spread research 
[4]. The concept of a quasi- invariant has been introduced in [5], and has been 
proven useful for representing the properties of nonlinear betatron dynamics. The 
quantum version of this invariant is developed, and from it quantum corrections to 
the classical results are found. The similarity between the Lie algebras associated 
with the classical case and the quantum case are used to obtain the relevant 
results. At first the classical linear and classical nonlinear cases are studied from 
the view point of their appropriate Lie algebras. These systems are quantized, and 
the corresponding Lie algebras are used to determine properties of the quantized 
systems. 

It is well known that the Courant-Snyder [6] invariant is particularly useful 
for determining the phase space pattern for the transverse dynamics of a particle 
in a storage ring. Using the Hamiltonian for a time-dependent simple harmonic 
oscillator, one can obtain the relevant invariant. However, when there are non- 
linear contributions to the Hamiltonian little success has been achieved in finding 
invariants. In order to better understand the behavior of a particle beam, it 
is useful to find an approximate invariant, which is associated with a nonlinear 
time-dependent Hamiltonian. The method used to find the quasi-invariant for the 
nonlinear oscillator is first used in the context of classical dynamics, based on us- 
ing the Lie algebra associated with elements obtained from powers and products 
of the position and conjugate momentum coordinates. To illustrate the method, 
an example is given for the linear system, where the invariant is exact, and the 
relevant algebra is SU(1,1) [7]. The method used for the linear system can be 
easily generalized to study a nonlinear one dimensional system. The method has 
the advantage that the time-dependent coefficients of the approximate invariant 
are found as the solution of a system of linear first order differential equations. 

For classical one dimensional transverse dynamics, an approximate invariant 
associated with a time-dependent Hamiltonian containing a nonlinear sextupole 
term is found. Both the Hamiltonian and the approximate invariant can be rep- 
resented as linear sums of the elements of a Lie algebra. The invariant is ap- 
proximate in the sense that terms of order greater than three, resulting from the 
Poisson bracket of elements of the algebra, are neglected. This results in a closed 
Lie algebra. The method is extended to quantum operators, and a similar Lie 
algebra involving operator elements is found. This is used to obtain the quantum 
quasi- invariant. The relation between the classical result is established with the 
aid of coherent states associated with the linear time-dependent oscillator. 

The Hamiltonian associated with nonlinear betatron dynamics studied in this 
paper is a special case of the Hamiltonian for a particle of mass m with the one 
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dimensional conventional form 



'•2 '•2 

H{t) = ^ + mu^oKit)^ + S{t)f. 1.1 

The method used to obtain this equation from a covariant formulation of storage 
ring dynamic is found in [2] and [8] . This equation is put in dimensionless form 

H{t) = ^ + K{t)^ + S{t)f, 1.2 

with the transformation to the dimensionless position, momentum, and energy 
variables q/qo, p/po, and H/Eq. Here 

= 1, 2 = %^ PoQo = —ri'^ 1-3 



mEQ ttujOq eo 

with eo = a;o^. In addition 



S(t) = 1.4 



1.5 



Introducing the dimensionless parameter 

e - ^ - i2- 
Poqa Eq 

gives for the quantum bracket of the dimensionless operators q and p 

[q,p\ = i ^ = ie. 1.6 

Poqo 

This allow the results which depend upon the quantum bracket to be expressed in 
terms of e. The quantum results, corresponding to various orders of h, are found 
with e = 1, and associated classical results are found with e — > 0, corresponding 
to the limit h 0. For applications to betatron dynamics, it is conventional to 
use go = 1 and po = \p\, which is the magnitude of the three-momentum of a 
relativistic particle. 

II. THE COURANT-SNYDER INVARIANT AND SU(1,1) 

The time-dependent Hamiltonian for one dimensional transverse dynamics is 
written in terms of the position coordinate q and the conjugate momentum p as 

H{t)^^ + K{t)^. 2.1 
The invariant Courant-Snyder associated with this Hamiltonian is 

/o(t) = Pit)P^ + Mt)PQ + lW ^ 2.2 



3 



which satisfies the partial differential equation 

The Poisson bracket of phase space functions /(p, q) and g{p, q) is defined as 

The functions ct(t), and 7(t) satisfy the equations 

da{t) 

dt 
d[3{t) 

dt 
dl{t) 

dt 

where 



= K{t)l3{t) - -f{t) 

= -2a{t) 2.5 
= 2K{t)a{t), 



Both the Hamiltonian Eq. (2.1) and the invariant Eq. (2.2) may be ex- 
pressed in terms of the elements of the Lie algebra SU(1,1). If one introduces the 
coordinates 

q + ip 



a = 



a = 



V2 

q — ip 



2.7 



with Poisson bracket 
then the functions 



V2 ' 

{a,a*}=i, 2.8 



Ai = c?, A2 = a*^, As = a*a, 2.9 
satisfy the Lie algebra of SU(1,1). Namely, 

{^1,^2} = 4z^3 

{^1,^3} = 2^1 2.10 
{A2, As} = -2iA2. 

In terms of the elements of the algebra Eq. (2.10), the Hamiltonian and the 
invariant become 

H{t) = ai{t)Ai + a2{t)A2 + a3{t)A3 ^ 
Io{t) = pi{t)A^ + P2{t)A2 + Ps{t)As. 



Requiring Io{t) to be real gives the relations 



2.12 



When these are substituted into Eq. (2.3), one finds, using Eq. (2.10), the set of 
linear differential equations 



' dt * 

dt 



dt 



'2ias{t) -2iai{t) 
-2ia3{t) 2iailt) 
4ia-i_{t) -4iai{t) 




2.13 



The functions ai{t) and Pi{t) satisfy the relations 

m- 



ai{t) = a2{t) = 
a3{t) 



1 



and 



ia{t) — 



4 

K{t) + 1 



i{t)-m 
1 



2.14 



2 

7W-^(t) 
2 



2.15 



m + 7(0 



These relations can be used to show that the system of linear differential equations 
Eq. (2.13) is equivalent to the system Eq. (2.5). With initial values given for /3(t) 
and d(3{t)/dt, the system of equations Eq. (2.13) can be integrated numerically, 
using Eq. (2.14) and Eq. (2.15). 
III. THE NONLINEAR SEXTUPOLE SYSTEM 



Next the method described above is extended to the classical case when a 
nonlinear term is added to the linear system Hamiltonian. However, in this case 
an approximation is made to obtain a finite closed Lie algebra which contains 
seven elements. As an example, one considers the Hamiltonian Eq. (1.2) where 
S (t) is the strength of the sextupole term [9] . Now defining functions of a and a* 





= a\ 


A2 


= a*2 


A3 


= a*a, 


A^ 


= a' 


^5 




Ae 


= a\* 


A, 


*2 

= a a, 
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one finds, keeping terms of order less than four in a and a*, the closed Lie algebra 



{^1 
{^1 
{^1 

{A2 
{A2 
{A2 

{As 
{As 
{A4 
{A4 
{A, 
{Ae 



A2}=4iAs 
A4} 
Ae} 
As} 
A,} 
Ay} 
A4} 
Ae} 
A,} 
Ae} 
Ag} 
A7} 



0, 

2iA4, 

-2^2, 

0, 

—SiA^, 
-iAa, 
0, 
0, 
0, 
0. 



{Ai,As} = 2iAi, 

{^1,^5} = 6^7, 
{^1,^7} = 4^6, 

{^2,^4} = -QiAe, 
{A2,Aq} = -4iAr, 

{A3, A5} = 32^5, 
{^3,^7} =^7, 

{^4,^7} = 0, 
{A5,A7} = 0, 



The Hamiltonian Eq. (1.2) may be written in the form 



i=l 



where ai{t), i 



3 are given by Eq. (2.14), and 
aQ{t) = a^it) = 3Q!4(t). 



3.2 



3.3 



3.4 



One can now find an approximate time-invariant associated with the Hamil- 
tonian Eq. (3.3). This is assumed to be of the form 



m = Y,mA 



J' 
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3.6 



which contains terms up to third order in a and a*. Since I{t) must be real, one 
finds 

When this along with the Hamiltonian Eq. (3.3) is substituted into Eq. (2.3), 
one finds, using the algebra Eq. (3.2), the system of linear first order differential 
equations 

d(3{t) 



dt 



= M(t)/3(t), 



3.7 
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where 



P{t) = 



3.8 



and 

M{t) = 

( 2ia3{t) 


-4^2 (^) 


4ia7{t) 





-2ia3{t) 
4iai{t) 


-2ia7{t) 
—6ia4{t) 
—4:iaQ{t) 



-2iai{t) 
2ia2{t) 


—Sia^lt) 
3^0:5 (t) 









6ia2{t) 









—3ia3{t) 


—6iai{t) 







-2iai{t) 


4^2 (^) 








2ia2{t) 
-4iai{t) 

-iOL2,{t) I 
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In these expressions the aii) * are given in Eq. (2.14) and Eq. (3.4). The first 

three /3(t) ^ are given in Eq. (2.15). The remaining /3(t) ^ are found as solutions 
to a system of first order differential equations Eq. (3.7). Using Eq. (3.5), the 
quasi-invariant may be written in the form 



/ =~ h{t) + 2$R(/?4(t)^4 + /?6(t)^6). 



3.10 



The first term /o(t) is the function Eq. (2.2), which is an invariant for the linear 
system. The remaining term may be expressed in the form 



with 



3.11 



ci(t)g^ + C2(t)g^p + C3(t)g/ + C4(^)p^ 

V^Ci(t) =3?/?4(t)+3?/?6(t) 

V^C2(t) = -(3$>/54(t) + '^P&it)) 
V2c3{t) = -(3$R/?4(t) - ^Peit)) 

V2c4it) = QP4{t)-^/36{t). 

The functions Ci{t), i — 1 ^ 4, satisfy the following system of first order differential 
equations: 

ci{t) = K(t)c2it) + 3Sit)ait) 
C2{t) = -3ci(t) + 2K{t)c3{t) + 3S{t)P{t) 
csit) = -2c2{t) + 3K{t)c4{t) 

C4{t) = -C3(t), 



3.12 
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where dot denotes differentiation with respect to t. 
IV. THE QUANTUM LINEAR SYSTEM 



As a first approximation for finding the quantum limits associated with the 
Hamiltonian Eq. (1.2) for transverse betatron oscillations, one neglects the non- 
linear multipole contributions and considers for each transverse degree of freedom 
a time-dependent harmonic oscillator with Hamiltonian 

H(t) = ?!±|W£, 4.1 

where t [c = 1) represents arc length along an ideal storage ring orbit. The 
dynamical evolution of the conjugate quantum operators p = q and q is determined 
from the Heisenberg equations 

The Cour ant- Snyder invariant as a function of the quantum operators q and p 
takes the form 

2/o(t) = [{wp - wqf + [q/wf]. 4.3 



The invariance follows from 



along with the conditions 

w + K(t)w- ^ = 

4.5 

'^ + K{t)q = 0. 

Expressed in the usual Cour ant- Synder parameters, one finds for each transverse 
coordinate q 

2/o = lit)f + a{t){qp + pq) + 4.6 

with 

a{t) = —WW 

Pit) =w^ 4 ^ 

^ Pit) • 
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The quantum states for this system can be constructed with the aid of the 
squeezing operator [10] and [11] defined as 

S = e , 4.8 

with complex ^ = \$^\exp{i4>) and boson operators a and d^. The time-independent 
rationahzed Hamiltonian is 

^o = ^^ = ata+l, 4.9 

where the boson operators a and d^ are found from 

a + a) ^ a — a) , 

with the commutation relations 

[q^p\=i =^ [a,d^] = l. 4.11 

The Courant-Synder invariant Eq. (4.3) or Eq. (4.6) is found from the time- 
independent Hamiltonian using the squeezing operator Eq. (4.8) to write 

h{t) = SHoS^= (^^^^+0, 4.12 



where 



h{t) = Se'^aS"^ = ]^(^ + w - iv^ a + ]^{^ - w - iv^ a) , 4.13 



with 

cosh 1^1 = -■\/{l/w + wY + w"^ 



tan 6* = ^ 4.14 

tan(6' -h 0) = 



1 — 

The eigenstates of Io{t) satisfy the eigenvalue equation 

io{t)\n,t) ={n+^ \n,t) 

kt) = ^|0). 
9 



4.15 



The states |n, t) are not Schrodinger states, for they are not solutions of the time- 
dependent Schrodinger equation 

i^^\n,t)s = H{t)\n,t)s. 4.16 

However, the Schrodinger states are of the form [12] 

\n,t)s = e''''^^\n,t), 4.17 
where the phase, as shown in Appendix A, is 

MO = -(n+i)/'^. 4.18 

To evaluate the quantum correction to /o(^) cind to find the uncertainties 
associated with the operators q{t) and one must use the appropriate coher- 
ent state associated with H{t). This state is the time-dependent generalization 
of the coherent state [13] obtained from the eigenstates of the time-independent 
Hamiltonian Eq. (4.9). This is the nearest quantum state to the classical state 
of the simple harmonic oscillator. The coherent state for a time-dependent simple 
harmonic oscillator can be generated from the squeezed ground state as 

= D(/3)|0,t)„ 4.19 

where the displacement operator D(/?) is defined as 

D(/3)=e^^'(*)-^*^(*). 4.20 

Here /5 is a complex parameter, which is the eigenvalue of the operator h{t). This 
parameter is related to the classical value of the invariant /o(^) since 

t\h{t)\i3, t)s = + i/mi\p\ = io{t) + ^/m- 4.21 

This includes the quantum correction (l/2)(^/|p|). The variance of an operator q 
is defined as 

cy\q) = s{(^.t\{q-q?W)s, 4.22 
where the mean value of the operator q is 

q = s{(3Aq\(^.t)s- 4.23 

Time-independent Hamiltonian Eq. (4.9) results are found using the coherent 
state I a), defined for the complex parameter a as 

\a) = D(a)|a) = e^^^-'^^'^lO), 4.24 
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where the parameter a is related to the classical coordinates of position q and 
momentum p as in Eq. (2.7). They are 

a{q) = a{p) = 

V 2 4.25 

h 

a{p)a{q) = -, 

which yield the minimum value for the uncertainty product. For the Hamiltonian 
Eq. (4.9), the coherent state \a) represents the quantum state nearest to the 
classical state, h ^ 0, for which a{q), cr(p), and the uncertainty product are zero. 

One can now use the states Eq. (4.19) and the definition of the variance Eq. 
(4.22) to obtain results appropriate for the a particle collider. For the scaling 
transformations frequently used in betatron dynamics 

p 

H ^ 4.26 

\P\ 

where the three-momentum magnitude \p\ is 

\p\^- 4.27 
c 

with relativistic particle energy £^ one finds that the uncertainties and the uncer- 
tainty product, represented in terms of the Courant-Snyder parameters Eq. (4.7), 
are 



a 



p\ rhi{t) 



a{q)a 



PJ V 4.28 

p dq 

dq\ _ h^mW) 



dt J 2\p\ 



Writing the amplitude as qamp — a/ (eo/7r)/3(t) with the emittance eo = 27r/o(t), 
one finds the results 



cr{qamp) _ /e 



'1 



1/2 



Qamp V ^0 / 4 29 

/ ^ ^particle 

2^~ 2^ 2~' 
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where e^/Tr, the quantum emittance, represents half the resolution distance of a 
particle in the beam. With the approximations 



hc^2x 10"^^TeVm 

4.30 

£ ^ 2TeV, 



one finds 

e^/TT 5 X 10"^° m. 4.31 
For a typical proton collider with P{t) ~ 300 m and with Qamp ~ 3.5 mm, one finds 

en/n « 4 X 10~®m 

-fi 4.32 
cr{qamp) ~ 3.9 x 10 ""mm. 



Similarly, the angular uncertainty is 



a ( ^ ] « 1.3 X 10-^^rad. 4.33 



V. THE QUANTUM NONLINEAR SYSTEM 

The method described can be extended to the case when a nonlinear term is 
added to the quantum Hamiltonian. As an example, one considers the quantum 
operator Hamiltonian Eq. (1.2). Defining operator elements of a and , with 
[a, = e, as 



il 


= a' is = 




As 


= {a^a + dd^)/2 A4 = d^ 






= ie = 


(a^a^ + aa^a + a^a^)/3 


Ar 


= (a+^a + d^dd^ + aa^^)/3. 





5.1 



one finds, keeping terms of order less than four in d and d^ and first order in the 
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quantum parameter e, the closed approximate Lie algebra 



[^1 


M] 


= 4ei3 


[Ai,As] 


= 2eii, 




M 


= 0, 


[Ai,A,] 


= Geiy, 


[il 


M 


= 2ei4, 


[Ai,A7] 


= 4ei6, 


[M 


M 


= -2ei2, 


[^2,^4] = 


= -Gels, 


[M 


M 


= 0, 


[^2,^6] 


= -4ei 


[A2 


M 


= -2ei5 






[h 


M 


= —3eA4, 


A, is] 


= Seis, 


[is 


M 


= -eie, 


[i3,i7] 


= eiy, 




M 


= 0, 








ie] 


= 0, 


[i4,i7] 


= 0, 


[h 


M 


= 0, 


[^5,^7] 


= 0, 




M 


= 0. 







5.2 



These algebraic relations are the same as those associated with the Poisson bracket 
relations. The hermitian Hamiltonian operator and the hermitian quasi-invariant 
operator are 



i=l 



and 



5.3 



5.4 



m^j2f^^(t)Aj, 

where the ai{t) and Pj{t) are defined as before. The equation which must be 
satisfied by the hermitian quasi-invariant operator is 



di{t) di{t) 



dt 



dt 



+ i[H{t),i{t)]=0. 



5.5 



This equation along with the quantum algebra leads to the same set of differential 
equations Eq. (3.7) that appear in the classical case. 
VI. QUANTUM CORRECTIONS 



The quantum corrections to the quasi-invariant are obtained using first order 
perturbation theory. The Boson operators which occur in the linear invariant Eq. 
(4.12) can be written as 



f2{w) 

fiH 



a 
at 



6.1 



13 



The inverse transformation is 
where 

h{w) = - i — + w - iw] 
2 \ w I 

1/1 ^'^ 

/2(«^) = o ( w-iw 

and 

|/iHp-|/2HP = l. 6.4 
These scaled Boson operators satisfy the commutation relations 

[a,at] = 1 
[6,6t] = i. 

The operators q and "p become 



6.5 



^= (t(;/v^)(S + 6t) 



6.6 



The quantum corrections to the quasi-invariant are found from the operator 

i(t) = io(t) + ii(t) 6.7 

where from Eq. (4.12) 

h{t) = {l/2)[{wp - wqf + {q/w)% 6.8 

and 



h{t) ^ci{t)f + C2{t)^{fp + qpq + pf) + 

C3it)^ip'^q + pqp + qf) + c^{t)f- 

The classical values of the operators q and p are 

q=s {P,t\q\P,t)s 
p=s {P,t\p\P,t)s. 



6.9 



6.10 



The expectation value of the quasi-invariant operator is 

s{P,t\iit)\P,t)s =s {P,t\ioit)\P,t),+s {P,t\h{t)\P,t)s, 6.11 

14 



where s(/5, t\Io{t) t)s is given in Eq. (4.21). The correction to the hnear invariant 
is 

s{P,t\h{t)\P,t)s^h + hq,. 6.12 
The classical correction to the linear invariant is 

h = ci(t)g3 + C2{t)fp + C3{t)p^q + C4{t)p^, 6.13 

and the quantum correction is 



hqc = {iq/2)[/3{t)ci{t)+l{t)cs{t) - a{t)c2{t)] 

+ {p/2mt)c2{t) + 37(t)c4(t) - a{t)cs{tmn/\p\). 



6.14 



VII. RESULTS AND CONCLUSIONS 

The nonlinear time-dependent Hamiltonian for one dimensional transverse 
classical dynamics is written in terms of the position coordinate q and the conjugate 
momentum p in Eq. (1.2). For this Hamiltonian, the equation of motion is found 
from Hamilton's equations 

. dH{t) 



P = 



dp 

dH{t) 



7.1 



dq 
to be 

q + K{t)q + ZS{t)q^ = 7.2 
The classical approximate invariant associated with this Hamiltonian is 

_ /3(t)p^ + 2c^(t)M + 7(t)g^ 

2 7.3 
+ ci{t)q^ + C2{t)q^p + cz{t)qp^ + c^{t)p^ . 

The time-dependent functions a{t), /3{t), and j{t) are found from the Eq. (2.5) or 
Eq. (2.13), and the functions Ci{t) can be found from the differential equations Eq. 
(3.12). These system of equations are equivalent to the system of linear equations 
Eq. (3.7). 

Numerical results are given which confirm the analytical development in the 
previous sections. Periodic solutions for the functions Ci{t) allow the determination 
of these functions at a fixed point in a lattice with a sextupole nonlinearity. The 
values of the functions q and p are determined from nonlinear tracking for the first 
five circuits of the lattice. After the j*'^ turn, the quasi-invariant becomes 

= Io{j) + ci{j)g{l,j) + C2{j)g{2,j) + C3{j)g{3,j) + c^AJ), 7.4 
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with 



gi2J)=qij)Mj) 



7.5 



Prom the requirement that 



/(/c) - /(I) = 0, 7.6 
for A; = 2 ^ 5, one finds the system of hnear equations 

4 

A(fc)=^Q(j)A^7(z,/c), 7.7 

with 

A(/c) = -(/o(/c)-/o(l)) 
^g{i:k) = g{i,k) - g{i,l) 

For the numerical results, the FODO approximation is used to find the lattice 
function f3{t). This function is derived in Appendix B. The system Eq. (7.7) and 
Eq. (7.8) is solved numerically to find the coefficients Ci{t). The thin lens approx- 
imation is used where the lattice is made up of a single thin sextupole element 
and identical cells of length L. Each cell consists of a focusing and defocusing 
magnet separated by a bending drift magnet. The focal length of the focusing and 
defocusing magnets is /, and the phase advance per cell ji is found from 

sin{^i/2) = A 7.9 

The tune v is obtained from 

.= t!k, 7.10 
27r 

where Nc is the number of cells. The maximum value of j3{t) occurs when a{t) = 0, 
and P{t) = 1/7 (t), and it is found from 

The phase space plot of PmaxP cm and q cm for the classical quasi-invariant 
/ is shown in Figure. 1. The classical results for both the invariant for the 
linear system Iq cm and the quasi-invariant / cm for the nonlinear system are 
plotted in Figure 2. as a function of turn number. It is clearly seen here that 
the methods leading to the quasi-invariant produce a more stable quantity than 
Jo- For the example considered, the values A^c = 4, = 7r/2, L/2 = 8875 mm. 
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and V = 0.33666667 + NciJ,/{27r), with near resonance fractional tune contribu- 
tion, have been used. The initial values q = 0.3 cm, (3^axP = along with the 
sextupole strength 3Se = 0.1 x 10~^ cm~^ have been used. For integer j, the 
sextupole function is approximated by S{t) — (se/3)5(t — JTq), where Tq is the 
orbital period. For the present case, the values of the periodic functions Ci (0) are 

ci = -3.41219 X 10-^ cm-^ 

-0.91910 X 10"^ cm-2 
+0.99563 X 10"^ cm-2 
-1.11468 X 10"^ cm-^ 

(7.12) 

with (3„^ax = 38389.279 cm. 

It is clear from the Figure 2. that the quasi- invariant is nearly stable. It 
remains this way for increasingly larger number of turns. It oscillates with small 
amplitude and with period of 100 turns. The amplitude of the oscillation de- 
pends upon the strength of the sextupole nonlinearity, and the period results from 
the nearness of the fractional tune to the third integer resonance. Although, the 
present quasi-invariant, which includes terms in q and p through third order, be- 
comes increasingly unstable for large values of the sextupole strength or large 
initial values of the amplitude q, it is clear that the method can be extended to 
include arbitrarily higher order corrections which will improve the stability of the 
quasi-invariant. The quantum correction associated with the quasi-invariant can 
be found from Eq. (4.21) and Eq. (6.14), and for the numerical example being 
considered it takes the value 

Igc ~ (1/2 + q/2{PmaxCl + C3/Pmax))n/\p\ ^ 0.538n/\p\. 7.13 

Although very small for a hadron collider, it would be more significant for a low en- 
ergy nonlinear time-dependent oscillator of the type described by the Hamiltonian 
Eq. (1.2). 

In conclusion, it is seen that the Lie algebra methods used for both the classical 
and quantum quasi-invariant s provide a useful approximation for the invariant as- 
sociated with the time-dependent nonlinear oscillator. For applications to betatron 
dynamics, this method provides a complimentary method to the usual nonlinear- 
map tracking methods. In addition, the quantum states \/3,t)s of Eq. (4.19) pro- 
vides the connection between the quantum operator for the quasi-invariant and 
the classical result when these states are used to form matrix elements of the type 
used to obtain the quantum uncertainties Eq. (4.28) and the quantum correction 
Eq. (6.14). 

This work was performed in part while the author was University Scholar 
in Theoretical Physics at UCLA and partially supported by U. S. Department 
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Institiiiid Teicneolaiochta Bhaile Atha Cliath grant 9571. Computations have 
been done with the aid of REDUCE and the CERN Computer Library. 
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APPENDIX A: THE SCHRODINGER STATE PHASE 



The phase Eq. (4.18) can be found by first differentiating Eq. (4.17) with 
respect to t, and then using Eq. (4.16) to write the matrix element 

+ {n,t\ — \n,t) = -i{n,t\il{t)\n,t) 

= ]-{w'^ +K{t)w'^ + l/w'^){n + 1/2), A.l 

where Eq. (4.13) is used to express the Hamiltonian Eq. (4.1) as a function of 
h{t) and its adjoint. The matrix elements of this operator can be found from 



b"'"|n,t) ^ Vn + l\n+l,t) A.2 

and 

b |n, t) = v^l'^ — 1) ^) =^ {n,t\h^ = ^/n{n — l,t\. A.3 
Making the replacement n — > n — 1 in Eq. (A.2), one can derive the identity 

{n,t\-^\n-l,t) + {n,t\h^ — \n-l,t) 

d 

= y/n{n,t\ — \n,t), AA 

1 

= -(i{ww — w^) — 2w/w)h + i{ww — w'^)hK A.5 
(n,t|^|n,t) = {0^t\-^^\0^t)+i^{ww-w^). A.6 



where 



Therefore 



Choosing 

(0,(|||0,*> = ^i=p^>, A7 

one finds from Eq. (4.5), Eq. (A.l), Eq. (A.6), and Eq. (A. 7) the differential 
equation 

dt ^ 2' w^{ty 
which has the solution Eq. (4.18). 
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APPENDIX B: THE BETA LATTICE FUNCTION 



In this appendix, following the methods of [6], the lattice function f3{t) used 
in the numerical calculations is derived. It is found for a lattice made up of 
similar cells of the FODO (focusing, drift, defocusing, drift) form. Focusing and 
defocusing are achieved with thin lens quadrupolc magnets, and the drifts occur 
through bending dipole magnets of length L and strength Bq. The function P{t) 
has period 2L, and the function on the interval L < t < 2L is found from that on 
the interval < t < L using 



The beta functions for a lattice with phase advance n per cell are found from 
the (1, 2) component of the transfer matrix. The function (3{t) is found from 



where the focusing and defocusing matrices for lenses of focal length / are, respec- 
tively. 



f{t)L<t<2L — fC^L - t)o<t<L- 



B.l 



Pit) sin = (0(t)FO(L)DO(L - t))i2. 



B.2 




B.3 



The matrix for a drift of distance t is 




BA 



The resulting beta function for < t < L is 



2L / L 1 L , 




B.5 



where L = 2/sin(///2). 
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FIGURE CAPTIONS 



FIG. 1. The quasi-invariant in phase space using PmaxP cm and q cm. 
FIG. 2. The Courant-Snyder invariant, eo/27r = Iq cm, and the quasi-invariant e/2n = 
I cm, as a function of turn-number. 
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Figure 1: The quasi-invariant in phase space using 
^max V cm and q cm. 
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Figure 2: The Courant-Snyder invariant, eo/27r = Iq 
cm, and the quasi- invariant e/27r = / cm, as a function of 
turn-number. 



